Non-Gaussian two-mode squeezing and continuous variable entanglement of linearly 
and circularly polarized light beams interacting with cold atoms 
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We investigate how entangled coherent states and superpositions of low intensity coherent states 
of non-Gaussian nature can be generated via non-resonant interaction between either two linearly 
or circularly polarized field modes and an ensemble of X-like four-level atoms placed in an optical 
cavity. We compare our results to recent experimental observations and argue that the non-Gaussian 
structure of the field states may be present in those systems. 
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I. INTRODUCTION 



Quantum state entanglement is a recognized resource 
for achieving efficient quantum communication protocols. 
With the emergence of proposals for continuous variable 
(CV) quantum communication protocols [H, Q , such as 
quantum teleportation d, [3| and dense coding [1, d, 01 , 
there has been an increasing interest in the generation 
and manipulation of entanglement in the CV regime 
through a diversity of experimental setups. Recently, bi- 
partite CV entanglement was achieved in a most remark- 
able experiment, through the interaction of a coherent 
linearly polarized light beam with a cloud of cold atoms 
in a high finesse optical cavity [1] . The end product of the 
interaction is the generation of two entangled squeezed 
modes with orthogonal polarizations. The entanglement 
between those beams is demonstrated by checking the 
inseparability criterion for CV The inseparabil- 

ity criterion is strictly a necessary and sufficient condi- 
tion for entanglement only for Gaussian states, but is a 
sufficient condition for entanglement for any other CV 
state. However it is not completely evident whether the 
bipartite state generated in this experiment is Gaussian 
or not. Precisely speaking the interaction between the 
two orthogonally polarized fields intermediated by the 
atomic cloud is highly expected to be nonlinear in the 
field modes annihilation and creation operators. As is 
well known only Hamiltonians which are at most bilinear 
in canonically conjugated variables can lead to Gaussian 
evolution [lli] . A deeper analysis of this system was given 
in a series of papers [l3l . [l4| . and some results were nu- 
merically confirmed [15| . However it would be certainly 
important to stress all the available possibilities for gen- 
erating entanglement and to infer on the Gaussian or 
non-Gaussian character in this experimental setup. 

In this paper we investigate a model for the inter- 
action of two orthogonally polarized quantum fields with 
an ensemble of X-like four level atoms, deriving an effec- 
tive Hamiltonian accounting for the field modes interac- 
tion. Up to first order the interaction results to be bilin- 
ear in the field operators, which however is multiplied by 
the difference of population of an effective ensemble of 



two-level atoms, possibly leading thus to a non-Gaussian 
evolution. Conditioned on the atomic population mea- 
surement, the two orthogonally polarized fields are left 
on an non-Gaussian entangled state, which shows similar 
properties of the Gaussian two-mode squeezed vacuum 
state. Since the experiment allows for two regimes de- 
pending on the detuning between the incident light and 
the atomic system, called self-rotation [l^ or polarization 
switching [8| , one deals either with a circularly polarized 
beam or a squeezed linearly polarized beam. By appropri- 
ately setting one of the input linearly polarized modes in 
the vacuum state, we obtain in the output a coherent su- 
perposition of coherent states in the orthogonally polar- 
ized linear modes. When viewed from the circular polar- 
ization frame, this superposition results in an entangled 
coherent state between the two modes in polarization + 
and — . The presence of this superposition, in one po- 
larization reference frame and an entangled state in the 
other can explain all the non-classical features observed 
in the before mentioned experiment, with a non-Gaussian 
state however. Remarkably recently much effort has been 
dedicated to the generation of such non-Gaussian states 
in propagating light fields by photon-subtraction [3 • 
In this paper we show that in principle such states may 
be generated in the experiments reported in 
as well. 
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This paper is organized as follows. In Sec. II we review 
the description of the interaction between the two orthog- 
onally polarized fields with an ensemble of N cold atoms 
following [Tsj and obtain an approximate solution to the 
Heisenberg-Langevin equations that govern the atomic 
ensemble evolution in the dispersive regime considered. 
We then derive an effective Hamiltonian accounting for 
the interaction of an ensemble of two-level atoms and 
fields in both circular and linear polarizations. We inves- 
tigate, in Sec. Ill, the dynamical generation of entangled 
coherent states, and conditional generation of superpo- 
sition states. In Sec. IV we analyze the squeezing of 
quadratures variances for both linear and circular polar- 
izations and infer on the inseparability of the two orthog- 
onal modes. Finally, Section V contains a summary and 
conclusions. 
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frequency is uj, the detuning from the atomic ensemble 
transitions resonance is equal to A = uJat — ^- The cou- 
pling constant between the atoms and field is g ~ eod/h, 
where d is the atomic dipolc and eg ~ ^hwjleV ^ where 
V is the volume of the cavity. The dipole decay rate 7, 
as depicted in fig. [TJ is decomposed in two orthogonal 
rates as 7 = 7|| + 7_l. Thus, the atoms-field interaction 
Hamiltonian is described by [20j 



where 



H = Ho + Hi 



Hq — EiCTii + i?2CT22 + £'30'33 + E^aa 



(2) 



(3) 



FIG. 1: (Color online) Scheme of the four level X-like atomic 
model for the 6S'i/2,-F = 4 to 6P3/2,F — 5 transition of 
Cesium. 



II. THE MODEL 

Light fields interacting with cold atoms can show a di- 
versity of interesting phenomena such as squeezing, two- 
mode squeezing and CV entanglement [J [l^, [3, [T^ . 
ft can also be employed for implementation of quantum 
logic operations [1^. Recent experiments have demon- 
strated the direct relation between one mode squeezing 
and two-mode entanglement under a linear (Bogoliubov) 
transformation of the fields polarization reference frame 
[1, 

[1^ . In those experiments, an a;- linearly polarized 
probe field is let to interact with a cloud of cold cesium 
atoms in a high finesse optical cavity. The cavity out- 
put x-polarized signal and the y-polarized vacuum are 
squeezed, which results in the entanglement of two or- 
thogonal circular polarization fields. The probe light field 
is red-detuned by about 50 MHz of the 651/2, F = 4 to 



and 



6-P3/27-P' — 5 transition. As previously argued this 



complicated transition can be modeled by an X-like four- 
level atomic structure. 

Following Ref. we consider the atomic system as 
being a set of N X-like four-level cold atoms in an opti- 
cal cavity driven by a linearly polarized field, as shown 
schematically in figure [1] We employ collective oper- 
ators to describe the N atoms ensemble (e.g. (T14 = 
X^i^i e*"*|l)j;(4|i), and we denote the operators related 
to circularly polarized fields by indexes a4_(_), which are 
defined from the standard linear polarization components 



and 



(1) 



We shall consider both orthogonally polarized linear 
modes Ux and Uy initially in coherent states and latter we 
assume the mode ay in its vacuum state. The atomic tran- 
sition frequencies are chosen in resonance {uat = W13 = 
LU24:) for simplification matters. In that case, if the field 



Hi = hg 



e-^"*a_a32 + e'"*a^_a23 



(4) 



The atomic evolution is appropriately governed by a 
set of quantum Heisenberg-Langevin equations , here 
given in a rotating frame with the probe frequency u) as 

(714 = -(7 + iA)cri4 - iga+(<Tii - 0-44) + Fi4, (5) 

(723 = -{1 + i^)(^Ti- i9a-{(T22- (^zz) + F2i, (6) 

(711 = 27_LCT33 + 2711 (744 - i.g(a+CTi4 - /i.e.) 4- Fii, (7) 

(722 = 27||(733 -f 27_L (744 - ig(aL(723 - /i.e.) + F22, (8) 

(733 = -27(733 -t- ig{a}_(T2z - h.c.) + _F33, (9) 

(744 = -27(744 + ig(a^(7i4 - /i.e.) + F44, (10) 

( with h ~ \) where Fij are the Langevin operators. 

We want to derive an effective Hamiltonian accounting 
for the interaction between two- orthogonally polarized 
field modes. For that we make the following assumptions: 

(i) We shall neglect the fluctuation on deriving a sta- 
tionary solutions to the above set of equations. That 
means that we shall neglect the Langevin operators by 
taking Fij = 0. 

(ii) To simplify the four level system of atoms to an 
effective system of two level atoms, we consider that 
A is very large (A ^ 7 ^ 5) in such a way that 
the higher levels 3 and 4 are not significantly populated 
(t733 = (744 = 0). Thus by taking Eqs. (|5I6|) in the sta- 
tionary regime ((714 = 023 = 0) results in (0) lowest order 
approximation 



(0) _ -iga+{t)uii{t) 



'14 



(7 + «A) 



and 



(0) _ -iga-{t)(T22{t) 



'23 



(7 + «A) 



(11) 



(12) 
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These solutions can be replaced back into Eqs. ((Tl fTU)) for 
the levels population. Rewriting the interaction term in 
a symmetrized form we obtain 



0-44 



27g^ I t , 1 /ON 



^27 



l+ (^2 + A2) (4«+ + l/2) 



fas 



(:^2^A^(«-«- + 1/2)^22 



1 



-27 

= 27j_cr44 



5 



(72 + A2) 

275' 



(ala_ + 1/2) 



^744, (13) 



c^sa, (14) 



(^2 + A2) + ^/^^^^^ 

^2 



+27 

O'll = 27j^cr33 



— + r 2^ A2n (°-°- + 1/2) 
7_L (7 + ^ ) 



^33,(15) 



-27 



21 

7J- 



^-2-^(a>+ + l/2)an 



0-44,(16) 



where we have already neglected the Langevin terms. In 
the stationary regime where the states 3 and 4 will be 
considered practically not populated we find 



f33 



CT44 



r 



(72 + A2) 



1 + 



5 



(72 + A2) 



(ala_ + 1/2) 



n -1 



x(a^La_ + l/2)cr22, 



(7^ 



(7^ 



-(a\aj, + 1/2) 



x(a^a+ + l/2)(Tii. 



(17) 



(18) 



Since we have assumed A 3> 7 3> 5, we must have 

2 

^ 1' ^'^^ from now on we shall only keep in 
([T7]) and (HH]) first order terms in this quantity. That is, 
we assume 



0'33 



CT44 



(72 + A2) 



(72 + A2) 



)(ala_ + l/2)a22, (19) 



{a\a+ + l/2)aii. 



(20) 



These stationary solutions for the populations of levels 
3 and 4 can be then included in the correspondent first 
order terms for the coherences between sates 1,4 and 2,3 
are as follows 



(1) _ -iga+jt) (c^ii - 0-44) 



'14 



and 



(7 + «A) 
(1) _ -iga-it) ((722 - 0-33) 



'23 



(7 + «A) 



(21) 



(22) 



We now are in position to consider the Heisenberg 
equations for the field operators 



igcri4, 

«.90-23, 



(23) 
(24) 



where cr 



14(23) — e 

solutions (PT]) and 



14(23)- Substituting the stationary 
we obtain 



g 



2(7 + iA) 



1 - 



2 t 
(72 + A2) 



g 



2(7 + iA) 



2 t 
5 a+a+ 



(7^ 



(25) 



a_ = — < luj - 



2(7 + iA) 



1 - 



2 t 



(7^ 



2(7 + iA) 



2 t 



(7^ 



(26) 



where we have defined cr^ = cm — cttii and we assumed 
(Til + "'22 ~ 1- Remark the presence of the non-linear 
terms proportional to a_|_a^_a+ and to a-a}_a-. They 
will lead to non-Gaussian evolutions over the field modes, 
i.e., if the initial field mode states are Gaussian they will 
be driven to non-Gaussian states. The same is true for 
higher order terms, which imply that invariably the field 
operators evolution will be non-Gaussian as inferred in 
the introduction. These nonlinear terms are at least pro- 
portional to g^; which from our assumptions is very small 
compared to the linear terms in g^. which alone may or 
not lead to Gaussian evolutions depending on the pre- 
pared atomic system state. Keeping only the terms up 
to linear in g^, Eqs. ([^5]) and (|26p simplify to 



g 



2(7 -hzA) 

g' 



g 



2(7 + iA) 



2(7 -hiA) 



g 



2(7 + iA) 



a+cT;, (27) 



a_cr„(28) 



which, in terms of the linear polarization operators are 
given by 



g 



2(7 -hiA) 

g" 

2(7 + iA) 



2(7 -hiA) 

r.2 



yrr, ,(29) 



2(7 + iA) 



ttxtTz. (30) 



Assuming a,^ 
these equations as 



x,y we can rewrite 



ig- 



2(7 + zA)''^'^^ ^'''^ ^^^~"2(7 + zA) 



a^(T^.(31) 



The appropriate effective interaction Hamiltonians cor- 
responding to the evolution of the circularly polarized 
modes corresponding to Eqs. (|27I28[) is 



H 



e// 



-iA(a'|_a 



a^a 



(32) 
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with 



where 



A = 



2(7 + zA) 



= Ai + z A2, 



Ai 



9^1 



A2 = - 



2(72 + A2)' 



2(72 + A2)- 



(33) 

(34) 
(35) 



Similarly for the linear polarization with respect to 
Eqs. (PT|) . we obtain 



H. 



{I) 



A (a], a.x 



a' a„ cr. 



(36) 



which is a bilinear coupling between the x and y polar- 
ization modes. Had wc considered the nonlinear terms as 
discussed above we would end up with nonlinear interac- 
tion term as well, such as aJ^axOyOj/, which would drive 
the system state to a non-Gaussian entangled state such 
as discussed in [2l|. 

Remark that, being A a complex number, both effective 
Hamiltonians are non-Hermitian. The non-Hermiticity 
comes from the spontaneous emission from states 3 and 
4 and from the fact that we have eliminated those states 
in our treatment, leading to non- unitary process. Our 
choice for the approximations (A 3> 7 3> (?) are in accor- 
dance with experimental values, which settle 7 « 2.6— 16 
MHz, A sa 130 - 327 MHz, and g around 2 Hz. To keep 
such relations between the parameters (A, 7, g), the ap- 
propriate relation for the real and imaginary part of the 
coupling of the effective Hamiltonians is IA2I ^ |Ai|. 
Thus the effective Hamiltonians assume the following 
form 



H. 



eff 



A2(a^^a_|_ - ala_)crj. 



and 



H^'jf = zA2(ata, 



4ay)(T, 



(37) 



(38) 



III. ENTANGLED COHERENT STATES AND 
SUPERPOSITIONS OF COHERENT STATES 



We shall take firstly the Hamiltonians (|32p and (|36p in 
full form and then we consider the limit |Ai| ^ IA2I to 
properly set the time scale. As we shall see shortly, these 
Hamiltonians will generate non-Gaussian states only in 
favorable situations for the population imbalance of the 
atomic ensemble. To take into account the possible tran- 
sitions between the two fundamental collective atomic 
states, we consider one of the two possibilities: 

(a) The ensemble of atoms is found in a coherent 
macroscopic superposition of states 1 and 2, \Xa) = 
(1/V2)(|l)at ± \2)at), where |l(2)),t = |11..1(22..2)). 



(b) Each atom of the ensemble is found in a coherent 
superposition of states 1 and 2: \x a) = ^iLi^{\^)i + 
|2).). 

Ensembles of atoms prepared in coherent superposition 
of states interacting with light are important for several 
applications in quantum optics and quantum communi- 
cation m, m, [23, m, m, 123, m and constitute an es- 
sential resource for generation of entanglement for the 
kind of light field states we analyze. Although the first 
situation (a) is not experimentally easy to achieve, it will 
illustrate more easily the final non-Gaussian nature of the 
field states. The second situation (b) is possible to be re- 
alized by applying non resonant classical light 7r-pulses to 
the atomic cloud, initially prepared in one of its ground 
states, before the interaction with the quantum field po- 
larization we are considering. Alternatively one could 
employ more recent techniques to generate the coherent 
superposition, such as the one proposed in Ref. [2^. The 
field state resulting in this case is a generalization of the 
previous one, also clearly non-Gaussian. 

(i) Circular Polarization. We consider both circu- 
larly polarized modes initially prepared in coherent states 
(|a)+ and |/3)-), where the subscript + or — designate 
the two orthogonal circular polarizations. Firstly we as- 
sume the atomic ensemble state (a) After some calcula- 
tion the time evolved atoms-field state is 



1 



m))c = ^[|ae-) + |/3e— )_|1),, 



zblae' 



-At\ 



+ l/3e^*)~|2)at] 



(39) 



which is obviously in an entangled state. Remark that 
the bilinear operation itself do not allow the two fields to 
be entangled if for example the atomic ensemble is found 
in one of the two states, \l)at or \2)at, or a mixture of 
them. The atomic superposition of states is an essential 
resource to generate entanglement. 

By conditioning the measurement of the atomic system 
in the same initial superposition (a), we find 



Tr+,.{\{ximm 

1 



|ae^*) + |/3e" 



At 



)_±|ae-^*) + |/30- 



(40) 



where 



N+ = 2 n ± e" 



cosh 2Ait (|Q|2 + |/3|^) cos 2A2t 



X COS a 



I sin2A2t]} 



(41) 



Since A has both real and imaginary parts we see that as 
the time evolves the two modes become more and more 
entangled, due both to an increase of the coherent states 
amplitudes and to a time dependent dephasing between 
the two components of the superposition. To illustrate 



5 



we set the time scale for |A2|^ = 7r/2, obtaining 
1 



l^('r/2|A2|)) 



± 




(42) 



Taking the A1/IA2I ~^ limit in Eq. \(p{tt/2\X2\))1 
z«)) + |H/?))_± I H«)) + 1 (43) 



assumes the following form 
1 



which is an entangled, but non-Gaussian state. 

To show that such non-Gaussian entangled state would 
be present in many experimental configurations let us 
consider the much simpler to achieve state (b). Follow- 
ing the same procedure, but for the second initial atomic 
ensemble state (b), we obtain after the conditioned mea- 
surement of the atomic system in the same initial ensem- 
ble of superposition states, the following state 



ixaim) 



N 



1 



(44) 



where Cf = jwkv- 



is the binomial coefficient, and 
N (t) the corresponding normalization factor. This is 



a two mode non-Gaussian entangled superposition state 
with + 1 terms in the superposition. Despite be- 
ing simpler to be experimentally realized this last state 
has a cumbersome structure than the previous case Eq. 

For illustration we want to keep the simplest non- 
Gaussian state, and thus from now on we consider only 
the results for the coherent superposition state (a), but 
the results should follow in similar fashion for the exper- 
imentally more accessible situation (b). 

(a) Linear Polarization. Under the same initial condi- 
tions (both field modes in the circular polarization pre- 
pared in coherent states and the atoms in state (a), we 
consider the effect of the evolution operator 



= exp [-a(ata, - ay)cr,t] . 



(45) 



for the linear polarization. This is a beam-splitter oper- 
ation with a phase conditioning on the atomic state. It 
is well known [s^, [Ml, Ull that the bilinear beam-splitter 
operation do not entangle classical states. As a conse- 
quence, if the states of mode (x) and (y) are prepared 
in coherent states they evolve as coherent non-entangled 
states, unless the atomic system is prepared in a super- 
position of the two ground states, 1 and 2. 

The two polarization frames are related by the field 
operators from Eq. (1). By employing this relation it is 
immediate that the two initial states are related by 

|a) + |/3>_ ^ \a')y\/3%, (46) 
where a' = i{a - P)/V2 and P' = {a + /3)/V2. Under 
the evolution (H5l) the above state takes the form 



\m)± = -^U^Ht)Dy{a')D,m\0)y\0)A\^)at±\2)at) 

= -^Dy (cosh (Ai)a' - i sinh (At)/?') (cosh {\t)(3' + i sinh {\t)a') |0)j,|0):,|l)at 
v2 

±-^Dy (cosh (At)a' + i sinh {\t)(3') (cosh {\t)[3' - i sinh (At)a') |0)j^|0):,|2)at 
v2 

= -^|cosh(At)a'-isinh(At)/3%(g) |cosh(At)/3'-t-isinh(At)Q;% |l)Qt 
v2 

±^|cosh(At)Q!' -Kisinh(Ai)/3'),y ® | cosh (At)/?' - i sinh (At)a').^ |2)at, (47) 
v2 

I 

where Di{a) is the displacement operator in a respective detection into the same superposition state (a) prepared 
to the polarization i = x,y. Conditioning the atomic at the beginning, we have 



I cosh {\t)a' - i sinh {Xt)(3')y ® \ cosh {\t)(3' + i sinh {\t)a')x 



^Tr^,-{\{ximW} ^ 

±1 cosh {\t)a' + i sinh {\t)l3')y ® \ cosh (At)/?' - i sinh (At)a') 
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Again for the choice \X2\t — 7r/2, we obtain 

1 



l^(V2|A2|))', 



± 



±- 



/7V± 
1 



«sinh(|^)a' 



cosh(gi)/3' 



(g) ismh(— — )/3' + cosh(— — )a' 
2A2 2A2 



(48) 



which for A1/IA2I turns out to be 



|^(7r/2|A2|))'± = 




(a + /3) 



V2 



.(a-J3) 



± 



' V2 



(49) 



To match further with the experimental conditions we 
set the a;-polarized mode in a vacuum state, such that 
a — /3 = 0, obtaining 



l^(V2|A2|))^ 



1 



V2a) ± 

V 



-V2a^^ ) |0),. 



(50) 

This state represents a superposition of coherent states 
for the y-polarized mode, namely the odd or even coher- 
ent state [33 for the + or - choice for the atomic state. 
This state is of a remarkable importance for many appli- 
cations on quantum information and computation [34| . 
Now, in the circular polarization Eq. (|43)) . the state fol- 
lows 

^ [|(^a)) + |(-^«))_ 



l'^(V2|A2|))^ 



±|(-*a)) + |(za))_] 



(51) 



and unlike the linear polarization case (superposition in 
one mode and vacuum in the orthogonal mode), the 
two orthogonal circular polarization modes are entan- 
gled. This is a particular realization of quasi-Bell states, 
the so-called entangled coherent states, as fully discussed 
in [3^, and employed in [2l[. 

To end this section we remark that through our calcu- 
lations we have neglected dissipative effects over the two 
field modes. Obviously dissipative effects will be always 
present and will drive the system to mixed states, reduc- 
ing the amount of entanglement. However, the inclusion 
of dissipative effects would not alter the non-Gaussian 
nature of the field modes. Thus, since we want to keep 
the evolution as close as possible from a Gaussian one, 
and to simplify our discussion, we will not consider the 
effects of dissipation for the field modes. 



IV. QUADRATURE VARIANCES AND 
ENTANGLEMENT CRITERIA 



Entanglement in the circular polarization can be in- 
ferred from non-classicality signatures on the linear po- 



larization (see e.g. j33. [35jV This is commonly realized 
by analyzing squeezing of the quadratures variances, as in 
Ref. |13| . Employing the usual definition of the quadra- 
ture operators 



X, 



-4) 



Y, 



1 

2^' 



(52) 



where I = +,—,x, or y for each of the circular or linear 
polarization mode, the variance for the Xi quadrature is 
given by (AXi)^ = {Xf) - {Xif, and similarly for Yi. 
These two quadratures in a polarization frame can be 
combined to indicate if there is entanglement in another 
polarization frame through the CV inseparability crite- 
rion [lo| . In terms of variance operators it is given 
by 

la^b = i[A2(X„ + Xb) + A''{Ya - Yb)] < 1. (53) 

For Gaussian states, Ia,b < 1 is a necessary and sufficient 
condition for entanglement. In our case where the states 
are non-Gaussian, while la.t < 1 undoubtedly indicates 
entanglement, when > 1 nothing can be said about 
the presence of entanglement. In this section we analyze 
the behavior of the quadratures variances as indicators 
of entanglement and the inseparability criteria for the 
non-Gaussian states obtained. These results should be 
compared to the ones of Fig. 2 from Ref. [l^, attributed 
to Gaussian states. 

(i) Quadrature variances for the circularly polarized 
modes. The quadrature variances calculated from the 
state of Eq. (UOl) , and with parameters fixed to reproduce 
the experimental situation |12l | in the case the x-polarized 
mode is in a vacuum state (a = f3) are depicted in Fig 
2. In FiglDJa) we see an oscillatory (periodic) behavior 
for both quadratures variances of the + circularly polar- 
ized mode. Notice that the variance for the quadrature 
yi(- is periodically compressed, oscillating bellow the ref- 
erence line at the value 1/4, while the variance for X^ 
always oscillates above this line. In Fig[2{b), the only 
modification relative to the previous case is the value of 
the ratio between the A's (A1/IA2I = 0.1) In that case the 
variance of the quadrature X^ is compressed as well with 
time. We see a similar behavior, but for longer times the 
amplitude of the oscillations increases due to the non- 
Hcrmitian nature of the Hamiltonian. The variances X_ 
and y_ show exactly the same behavior of the ones for 
the -|- polarization. 
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FIG. 2: Time evolution of the quadrature variances for 
circularly polarized modes, with /3 = a, _Re[Q] — and 
Im[a] = 0.3. (a) (AX+)^ (continuous line) and (AY+f 
(dashed line) for A1/IA2I = 0. (b) (AX+)^ (continuous line) 
and (Ay+)2 (dashed line) for Ai/|A2j = 0.1. 

For the linearly polarized modes the criteria writes 

+A^iY,-Yy)ia,P)] < 1, (54) 

and for the situation a ~ P the inequahty correctly 
indicates no entanglement, since the x mode is left in 
a vacuum state. The plots for this case are shown in 
FigH FigUa) is for A1/IA2I = 0, while fig. 3(b) is for 
•^i/|-^2| = 0.1. In all situations the inequality is violated, 
which in this case correctly indicates no entanglement. 

(ii) Quadrature variances for the linearly polarized 
modes. For the linearly polarized modes we shall consider 
the variances from the state of Eq. PO]) with a = f3, and 
A1/IA2I — * 0. As can be noted in the plots of FigHfa) for 
this case and A1/IA2I = 0, only the variance of the Xy 
quadrature is squeezed (always below the reference line), 
while the x-mode is left in a vacuum state. Again in Fig. 
4(b) we take A1/IA2I = 0.1 and remark that the variance 
of the quadrature Yy is squeezed as well with time as an 
effect of the non-Hermiticity. 

Since the x polarized mode is always in a vacuum state 
the variances of its quadratures will not change with time. 
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FIG. 3: Plot of the time evolution of the inseparability cri- 
teria for the linear polarization for q = /3 fixing three dif- 
ferent lines for 7?e[a] = 0, Im[a] — 0.3 (continuous line), 
Im[a] — 0.7 (dashed line), and /m[Q] = 1.5 (short-dashed 
line), (a) A1/IA2I = and (b) Ai/iA2| = 0.1. All plots indi- 
cate no entanglement in accordance with the separable, but 
not Gaussian, state obtained in this situation. 

Thus the fact that one of the variances of one of the 
modes is squeezed bellow the noise limit is a good indi- 
cator that entanglement may be occurring in the circular 
polarization. Indeed if we plot the inseparability criterion 
for the circular polarization 

+A2(r+-r_)(a,/?)], (55) 

where the variances are obtained using the previous ex- 
pressions for this polarization, we observe a clear signa- 
ture of entanglement, as can be seen in FigO We remark 
however that there are situations (parameter choices) 
where the criteria will not indicate entanglement, while 
the state is clearly entangled. This is not a surprising 
fact since the state is non-Gaussian and thus violation of 
criteria is not a necessary ingredient for the existence of 
entanglement. For comparison we plot in FiglHlthe one 
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FIG. 4: Time evolution of quadrature variances for linearly 
polarized modes, with a = P, Re[a] = and /m[a] = 0.3. (a) 
{AXyf (solid line) and {AYyf (dashed line) for A1/IA2I = 0. 
(b) (AXyf (solid line) and (AYy)^ (dashed line) for A1/IA2I = 
0.1. Mode X is left in a vacuum state. 



mode reduced linear entropy, S{a, f3) = 1 — Tr{p^} time 
evolution, which in this case, since the joint system is 
pure, indicates entanglement between the modes. It is 
clear to notice that the criterion is correctly indicating 
entanglement for the same situations considered in Fig. 

m 

It is interesting to remark that the points of Figs. 
and [5] that indicate maximal entanglement are given by 
the state (|5ip for the circular polarization, which on its 
turn corresponds to a superposition of small amplitude 
coherent state ([5(7)) for the linear polarization. As we 
discussed this state has presented squeezing of one of its 
quadratures variance. Now we want to question how far 
is this state from the vacuum squeezed state, commonly 
attributed as being the case in many experimental situa- 
tions? As a matter of fact, by following the discussion in 
[13, the fidelity, F(|0,Pcat) = \/{^\Pcat\0, between this 
superposition state in the y polarization and the squeezed 
vacuum state. 



00 



-ie'^anhdei) 



\2n), (56) 
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FIG. 5: Plot of the inseparability criteria as a function 
of time for the circular polarization modes for a = (3 fix- 
ing three different lines for Re[a\ = 0, Jm[a] = 0.3 (solid 
line), Im[a\ = 0.7 (dashed line), and Im[a\ = 1.5 (short- 
dashed line), (a) A1/IA2I = 0. Curve points bellow 1 clearly 
signal entanglement of the circularly polarized modes (b) 
A1/IA2I = 0.1. The criterion is affected by the non-Hermiticity 
of the Hamiltonian, but at initial times it is still robust sig- 
naling entanglement, as depicted in the inset. 



can be very close to 1 as depicted in Fig. [TJ and in some 



sense they are "similar" states. For small values of a 
(which represents the states generated in the mentioned 
experiment), the fidelity (with the squeezed vacuum state 
with small squeezing parameter) is very high. However, it 
is important to notice that the superposition state is not 
a Gaussian state, although the squeezed vacuum state is 
Gaussian. The fact of being or not Gaussian is important 
for the criterion of entanglement used to interpret the 
results in Josse et al [3] ■ Moreover the Entanglement of 
Formation [sj expression for symmetric Gaussian states, 
employed to quantify entanglement in Ref. p^ . in this 
case is only a lower bound for the entanglement of the 
two modes. 
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FIG. 6: Plot of linear entropy time evolution for the + cir- 
cular polarization mode for a = P fixing three different lines 
for _Re[a] = 0, Im[a\ = 0.3 (continuous line), /m[Q] = 0.7 
(dashed line), and /m[Q] = 1.5 (short-dashed line), (a) 
A1/IA2I = Indicates a periodic entanglement and disentan- 
glement of the two orthogonal modes, (b) A1/IA2I = 0.1 The 
periodicity of entanglement-disentanglement is affected by the 
non-Hermiticity of the Hamiltonian. The system tends to be 
highly entangled as the time goes on. 



V. CONCLUSION 

In this paper wc have investigated a model [T^ for 
the interaction of two quantum fields with an ensemble 
of X-like four level atoms, and we derived an effective 
Hamiltonian accounting for the field modes interaction. 
We have demonstrated that the Heisenberg equations for 
the field modes operators are non-linear leading to a non- 
Gaussian evolution. This non-linear evolution will lead 
to an entangled non-Gaussian state or to a non-Gaussian 
superposition of coherent states, when viewed from the 
circular or linear polarization reference frame, depending 
on the initial states for the two modes and atomic sys- 
tem. Even when the evolution is kept as close as possible 
from a Gaussian one, i.e., bilinear in the two field modes 



FIG. 7: (Color online) 3D plot of the fidelity between the generated 
state fEg lSOI I and the squeezed vacuum state (Eq[56} as a function 
of a and The coupling condition adopted for the plot is A1/IA2I ~ 
0. 



operators, a superposition state of two atomic degenerate 
fundamental collective states can lead to a non-Gaussian 
evolution. By appropriately setting one of the input lin- 
early polarized modes in the vacuum state, we obtain in 
the output a coherent superposition of coherent states in 
the orthogonally polarized linear modes. Although this 
state is non-Gaussian it preserves similarities with the 
one-mode squeezed vacuum state. When viewed from 
the circular polarization frame, this superposition results 
in an entangled coherent state between the two modes 
in polarization + and — . We have compared qualita- 
tively these results to recent experimental results 
with a similar system, which however have attributed a 
Gaussian nature to the Quantum fields. The presence 
of this superposition, in one polarization reference frame 
and an entangled state in the other can explain all the 
non-classical features observed in this experiment, with 
a non-Gaussian state however. 

Linearization procedures for field operators are com- 
mon in quantum optics whenever non-linear processes are 
present. While it can lead to a good approximation to 
non-classical features such as squeezing, it does not allow 
to correctly infer the non-Gaussian nature of the system. 
By linearizing one is at most mimicking the second order 
moments which are present in the most general nonlinear 
case. Thus it is commonly attributed a Gaussian evolu- 
tion to the system as well, which as we have showed it is 
not particularly correct. Inference on quantum properties 
of single systems may well be described by the lineariza- 
tion procedure, but the description of entanglement of 
(at least) two modes suffers from inconsistencies which 
are only resolved when dealing with the correct nonlin- 
earized case. One typical feature that cannot be cor- 
rectly inferred is the separability of the two modes, since 
the criterion employed is only necessary and sufficient for 
Gaussian states. We have shown that invariably in the 



10 



present system the quantum state is non Gaussian. Re- 
markably recently much effort has been dedicated to the 
generation of such non-Gaussian states in propagating 
light fields by photon-subtraction [13, [3 ■ '^'^ demon- 
strated, in principle such type of states may begener- 
ated in the experiments reported in [1, [13, [H, and 
any other similar situation, whenever the atoms can be 
prepared in coherent superpositions if nonlinearities are 
to be neglected. If the nonlinearities are present, the en- 



tangled two mode field state will always be non-Gaussian 
independently of the atomic system state preparation. 
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